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Abstract. We discuss the slow, nonequilibrium, dynamics of spin glasses in their 
glassy phase. We briefly review the present theoretical understanding of the spectacular 
phenomena observed in experiments and describe new numerical results obtained in 
the first large-scale simulation of the nonequilibrium dynamics of the three dimensional 
Heisenbcrg spin glass. 



1. Why do we study spin glass dynamics? 

Spin glasses can be seen as one of the paradigms for the statistical mechanics of 
impure materials. Experimentally, however, the spin glass phase is always probed via 
nonequilibrium dynamic experiments, because the equilibration time of macroscopic 
samples is infinite. Simulations can probe equilibrium behaviour for very moderate 
sizes only, so that the thermodynamic nature of the spin glass phase is still a matter of 
debate. It is also as a model system that the glassy dynamics of spin glasses has been 
studied very extensively in experiments, simulations, and theoretically in the last two 
decades 0121- Although many theories account for the simplest experimental results, 
such as the aging phenomenon, early experiments revealed several other spectacular 
phenomena (rejuvenation, memory, etc.) that are harder to explain, allowing one to 
discriminate between various approaches 

In recent years, several theoretical descriptions of the slow dynamics of spin glasses 
described the physics in terms of a distribution of length scales whose time, t, and 
temperature, T, evolution depends on the specific experimental protocol, as reviewed 
in Ref. Aging is described as the slow growth of a coherence length, irit), 

reflecting quasi-equilibrium/nonequilibrium at shorter/larger length scales. Sensitivity 
to perturbations of quasi-equilibrated length scales accounts then for rejuvenation 
effects, while the strong temperature dependence of the growth law irit) explains 
memory effects [H 03 EE IE] • If early numerical studies revealed the existence of such 
a distribution of length scales its physical relevance was critically discussed only 
relatively recently [101 HI] . A major problem, however, is that most studies focused on 
the Edwards- Anderson model of an Ising spin glass, defined by the Hamiltonian 

H = - J2 JijSiSj, (1) 

<i,j> 
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where Si = ±1, the sum is over pairs of nearest neighbours of the chosen lattice, and Jjj 
is a quenched random interaction, drawn from a symmetric distribution. The interest of 
the model (JTJ) is that it has been studied very extensively so that some — but only very 
few! — issues have been settled, most notably the existence, for space dimensions d > 3, 
of a second order phase transition to a spin glass phase [T2~| IT3] . The nonequilibrium 
dynamics of the Ising spin glass has also been quite extensively studied. Unfortunately, 
for d = 3, some of the key experimental observations are not reproduced ITlj . 
although simulations in d = 4 have been more successful [H]. This may not be too 
surprising, since real spin glasses are made not of Ising spins but vector spins. When 
the interaction between spins is isotropic, the system is therefore best described by the 
Heisenberg spin glass Hamiltonian 



where the Sj are now three-component vectors of unit length. The Heisenberg spin glass 
has been far less studied than the Ising one, both statically and dynamically, presumably 
because it was hoped that the understanding of the apparently simpler Ising case would 
be sufficient to interpret experiments. Very recent experiments systematically comparing 
Ising (i.e. very anisotropic) and Heisenberg samples have shown substantial quantitative 
differences between the two types of samples, the nonequilibrium effects being indeed 
much clearer in Heisenberg samples jT^l EI] ■ 

Hence, it can reasonably be hoped that dynamic studies of the Heisenberg spin 
glass in d = 3 will reproduce the key experimental effects, so that deeper theoretical 
knowledge of the nature of the nonequilibrium dynamics of spin glasses can be gained. 
In this paper, we extract some preliminary results from the first large-scale numerical 
simulation of the nonequilibrium dynamics of the three-dimensional Heisenberg spin 
glass [T7j . 

2. Simulation details 

We simulate the Heisenberg spin glass (j2J in d = 3. The sum in (j2j) runs over nearest 
neighbours of a cubic lattice with periodic boundary conditions. We use a heat-bath 
algorithm JH] i n which the updated spin has the correct Boltzmann distribution for 
the instantaneous local field. This method has the advantage that a change in the spin 
orientation is always made. We use a rather large simulation box of linear size L = 60, 
and study several temperatures T = 0.16, 0.15, 0.14, 0.12, 0.10, 0.08, 0.04 and 0.02. 
Although all the quantities we shall study are self-averaging, we use several realizations 
of the disorder, typically 15, to increase the statistics of our data. 

Contrary to the Ising case, the phase transition of the Heisenberg spin glass is 
still an open problem. A decoupling between spin and chiral degrees of freedom was 
theoretically suggested while early simulations even questioned the mere existence 
of a phase transition ^H]- Very recent simulations involving the most efficient tools 
used to study the Ising spin glass, conclude that the model is characterized by a phase 
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Figure 1. Time evolution of the energy density J3J after a quench from infinite 
temperature at the initial time t w for T = 0.16, 0.15, 0.14, 0.12, 0.10, 0.08, 0.04 and 
0.02 (from top to bottom). 



transition, at T c ~ 0.16, where both spins and chirality simultaneously freeze [20]. This 
motivates our choice for the upper temperature studied in our dynamical approach and 
our use of spin variables as dynamical objects of study. 

3. What has to be measured? 

Before embarking in the complex phenomenology of spin glasses, it is worth discussing 
the simplest protocol one can think of to probe the spin glass phase. A 'simple aging' 
experiment consists of a sudden quench at initial time t w — from a temperature well 
in the paramagnetic phase, T 3> T c , to a constant, low temperature below the spin glass 
transition, T < T c . Aging means a very slow evolution with time t w (called 'age') of 
the physical properties of the system. To study this behaviour, we record two types of 
quantities. First, 'one-time' quantities can be studied, such as the energy density of the 
spin glass, 



The time evolution of e(t w ) for various low temperatures is presented in Fig. ^ from 
which the slow decrease of the energy towards an asymptotic equilibrium value is indeed 
observed, the sign that the dynamics is non-stationary. 

We also study 'two-time' dynamic quantities. While experiments usually record 
response functions, it is easier to measure the corresponding correlation functions in 
numerical work. Here, we record the spin-spin autocorrelation function defined as 



The qualitative behaviour of this function is well-known, and a prototypical example is 
shown in Fig. 121 As usual, the time decay of C(t + t w , t w ) can be decomposed into two 
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Figure 2. Spin-spin autocorrelation function J3J as a function of the time difference 
t for various t w logarithmically spaced in the interval t w £ [2, 57797]; t w increases from 
left to right. The temperature is T = 0.14. 



parts. For short time separations, t <C t w , the dynamics is almost independent of t W} 
while the later decay, t ^> t w , becomes slower the larger t w . Non-stationarity is reflected 
in the fact that C(t + t w ,t w ) ^ C(t). The physical interpretation is simple: since the 
relaxation time of the sample is infinite, the only relevant time scale is the age of the 
sample t w which imposes an age-dependent relaxation time: the older the sample, the 
slower its relaxation becomes. A careful analysis of short and large time behaviours of 
C(t + t w ,t w ) and comparison to experimental data is described in Ref. |17j . 



4. Understanding aging in real space 

The key problem is to understand the subtle slow changes that the system undergoes: 
what does 'old' or 'young' really mean for the sample? The answer necessarily connects 
to equilibrium, since the system eventually equilibrates for t w —>■ oo. Moreover, the 
decomposition of the decay C(t + t w ,t w ) between a fast stationary process and a slow 
non-stationary one directly suggests the existence of some sort of local equilibrium 
within the sample: a spin appears locally equilibrated (short-time dynamics) although 
the sample as a whole is still far from equilibrium and evolves towards equilibrium 
(long-time dynamics). 

It is possible to illustrate this last statement, as was done in the Ising case 
Because of the disorder, the spin orientations in an equilibrium configuration are 
random, so that it is impossible to detect any domain growth by simply looking at 
the spin directions. However, two systems, (a, b), evolving independently but with the 
same realization of the disorder will reach correlated equilibrium configurations , so 
that the orientation of the spins in the two copies will be similar, up to a global rotation. 
In Fig. |3J we present pictures where the 'orientation' variable 

ooB^(t) = S?(t)-Sj(t) (5) 
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Figure 3. The orientation variable cos 8i denned in Eq. © is encoded on a grey scale 
in a 60 x 60 x 60 simulation box at three different times t w = 2, 27 and 57797 (from 
left to right) and temperature T = 0.04. The growth of a local random ordering is 
evident. 



is encoded on a grey scale. Comparing three successive times, it becomes clear that 
aging involves the growth with time of a local random ordering imposed by the disorder 
of the Hamiltonian. 

It is of course possible to go beyond simple pictures of black and white domains 
and measure the growing coherence length £r{t) corresponding to the mean domain size 
in Fig. 01 For this purpose, the spatial decay of the following correlation function is 
recorded: 

C 4 M) = -^£S?(f) ■ S- +r (t) S}(f) ■ St +r (t). (6) 

This function is a straightforward generalization of the two-spin, two-replica correlation 
function studied in the Ising case which measures spatial correlations of the random 
relative orientation of two spins P|. In Fig. 01 we show this function for the same 
parameters as the correlators of Fig.EJ The spatial decay of C 4 (r, t) is clearly slower for 
larger t, in agreement with the pictures in Fig. EJ Physically, this means that a larger 
time t w implies a slower relaxation due to a larger coherence length, very much as in 
standard coarsening phenomena. 

Note that due to periodic boundary conditions, the function (0) in Fig. HI is 
symmetric about L/2 = 30. In Ref. [TH], it was argued that spin and chirality degrees 
of freedom undergo different aging dynamics because they are statically decoupled. The 
numerical support for this statement was the observation, for a system of linear size 
L = 15, that the autocorrelation (jlj becomes stationary at large t w . From Fig. we 
immediately recognize that the data of Ref. [TH] are plagued by severe finite size effects, 
so that the conclusions of previous aging studies of the Heisenberg spin glass [T[JJ 121] 
must be taken with some care and justifies our numerical effort of simulating a very large 
system, L = 60. The scaling properties of C^(r,t) and the properties of the coherence 
length are further discussed in Ref. |17j . We make here the important remark that 
much larger length scales can be reached in the same numerical time window for the 




Figure 4. Two-spin, two-replica correlation function © as a function of the distance 
r between the spins for various t w logarithmically spaced in the interval t w £ [2, 57797]; 
t w increases from left to right. The temperature is T — 0.14. 



Heisenberg spin glass than for the Ising case, which may indicate that richer behaviour 
can be seen in non-equilibrium simulations of the Heisenberg spin glass than the Ising 
one. 

5. Conclusion 

We have motivated the need for large-scale numerical simulations of the three 
dimensional Heisenberg spin glass in order to fill the gap between spatial theoretical 
descriptions of spin glass dynamics and experimental observations. The results presented 
here for the dynamics of the model (J2J show that spin variables qualitatively follow the 
same type of aging behaviour as in the Ising case, which is due to the slow growth 
with time of a dynamic coherence length. In Ref. we analyze in detail the scaling 
properties of the dynamic functions reported here. The observation that very large 
length scales can be reached in the numerical time window, see Fig. 0J gives us the 
hope, also confirmed by preliminary work, that the model will allow us to reproduce most 
of the experimental effects, with the advantage that simulations have direct access to 
the distributions of length scales involved in phenomenological theories, gaining further 
understanding of spin glass dynamics. 
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